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We study the Okamoto-Nolen-Schiffer (ONS) anomaly in the binding energy of mirror nuclei at 
high density by adding a single neutron or proton to a quark gluon plasma. In this high-density 
limit we find an anomaly equal to two-thirds of the Coulomb exchange energy of a proton. This 
effect is dominated by quark electromagnetic interactions — rather than by the up-down quark mass 
difference. At normal density we calculate the Coulomb energy of neutron matter using a string-flip 
quark model. We find a nonzero Coulomb energy because of the neutron's charged constituents. 
This effect could make a significant contribution to the ONS anomaly. 
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The Okamoto-Nolen-Schiffer (ONS) anomaly is the 
long-standing discrepancy between the calculated and 
measured binding-energy differences of mirror nu- 
clei . The anomaly is likely to arise from charge sym- 
metry breaking (CSB) in the strong interaction ||, itself 
believed to originate from the up-down quark mass differ- 
ence and electromagnetic effects in the Standard Model. 
Thus, the study of CSB is a useful tool to elucidate the 
structure of strongly-interacting nuclear systems. 

The ONS anomaly can be calculated on several lev- 
els. Perhaps the simplest is the observation by B. A. 
Brown Q that the magnitude of the anomaly is approx- 
imately equal to the Coulomb exchange energy. If one 
adds an extra proton to a nucleus in a simple Hartree- 
Fock picture, there will be both a direct (Hartree) and 
exchange (Fock) Coulomb interaction with the other pro- 
tons. If one — arbitrarily — neglects the Fock term, one 
obtains a better agreement with experiment. 

At a different level, Blunden and Iqbal compute the 
ONS anomaly by calculating the contribution from p-to 
mixing to the CSB component of the nucleon-nucleon 
(NN) interaction ||. Their CSB interaction can explain 
part of the anomaly. However, at present this is con- 
troversial, both in the choice of meson couplings || and 
in the momentum dependence of p-u) mixing §. There 
have been also a number of calculations of the anomaly 
based on the contribution from the up-down quark mass 
difference (Am). Indeed, Nakamura and coworkers ||] 
calculate a CSB NN interaction using a constituent quark 
model where the short-range color hyperfine interaction 
depends explicitly on the quark masses. Moreover, the 
mass difference between the neutron and proton may 
be density dependent [[1. Finally, there are other more 



recent model calculations, such as the one reported in 
Ref. 0. 

Although the observation by Brown is not a dynamical 
explanation, it is an interesting characterization of the 
size of the anomaly. Could there be something wrong 
with the exchange term? As the nucleon is a composite 
object, could it be that the exchange energy of composite 
objects yields results significantly different from the ex- 
change of point nucleons? One expects identical results if 
the composite scale of the nucleon is much smaller than 
the inter-particle spacing. However these scales are sim- 
ilar in nuclei. Moreover, although various calculations 
based on the up-down quark mass difference exist, we 
are not aware of any calculation of electromagnetic (EM) 
effects between quarks to the ONS anomaly. 

The neutron-proton mass difference in free space is 
made up from comparable contributions of A771 and EM 
effects. Note that EM and Am terms contribute with op- 
posite signs to the neutron-proton mass difference. How- 
ever, the ONS anomaly is sensitive to the density depen- 
dence of these contributions so the relative sign is un- 
known. In this letter we study EM effects involving the 
Coulomb exchange interactions of composite nucleons. 

To clarify the importance of EM and Am terms we 
consider a high-density limit of the ONS anomaly. We 
will show, with some mild assumptions, that in the high- 
density limit: (1) there is an ONS anomaly and (2) that it 
is dominated by EM effects with Am being unimportant. 
Further, (3) the magnitude of the anomaly is simply re- 
lated to the Coulomb exchange energy and (4) its sign is 
the same as that observed at lower densities. Finally, we 
will perform model calculations to see how relevant this 
high-density limit is to normal-density nuclei. 
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Consider very high-density symmetric nuclear matter. 
We assume that an electron gas makes the system elec- 
trically neutral. Thus the direct Coulomb interaction 
vanishes. Yet Coulomb exchange effects are still present. 
Now add either one proton or one neutron to the system 
and calculate the change in energy. First, model the sys- 
tem as a Fermi gas of elementary nucleons. An added 
proton will have a Coulomb exchange energy of 

V p = -e 2 ^ . (1) 

TT 

Here kp is the Fermi momentum of the proton and e 
is its electric charge. In contrast, an added neutron has 
zero Coulomb exchange energy: V n = 0. Thus, the energy 
difference between an added proton and a neutron is just: 

E p — E n = V p + M p — M n = — e 2 — - — AM , (2) 

TT 

where AM = M n -M p = 1.29 MeV is the neutron-proton 
mass difference. Equation (g) is the simple expectation 
of a model with unexcited point nucleons. 

Next we consider a quark-gluon plasma. We assume 
because of asymptotic freedom, that at very high density 
the system is nearly a free Fermi gas of quarks. This 
is because the strong coupling ots(kp) becomes small at 
the large momentum scale characterized by kp . When 
a proton is added it will dissociate into two up and one 
down quark. Therefore, the Coulomb exchange energy of 
these three quarks is 

where denotes the quark electric charge and kp is 
the quark Fermi momentum. Note that there are three 
times as many (valence) quarks as nucleons. However 
the quarks have an extra color degeneracy of three. As 
a result, the quark Fermi momentum in Eq. (0) is the 
same as the proton's Fermi momentum in Eq. (pi). The 
sum of the squares of the valence charges in a proton is 
(4/9 + 4/9 + l/9)e 2 = e 2 . Because of this "numerical ac- 
cident" the quark Coulomb exchange energy is equal to 
the Coulomb exchange energy of an elementary proton. 

An interesting difference arises when we add a neutron. 
In a quark-gluon plasma the Coulomb exchange energy is 
no longer zero because a neutron is made up of charged 
constituents. Moreover, the exchange energy is always 
negative independent of the sign of the charges; the con- 
tributions from positive and negative charges add rather 
than cancel. Indeed, the sum of the squares of the valence 
quark charges in a neutron is (4/9+l/9+l/9)e 2 = 2e 2 /3. 
Thus, the neutron Coulomb energy is fully two thirds of 
that of a proton: Vn q ^ = —^e 2 kp/n. The energy differ- 
ence between an added proton and a neutron becomes: 

Bta) _ B M = -e 2 ^ + -e 2 ^ = --e 2 ^- . (4) 

F ' TT 3 TT 3 TT 



We choose to define an ONS anomaly Ai?oNS a s the 
actual energy difference, which we assume is given by 
Eq. (|4|), minus the hadronic- model expectation of Eq. (||) 

ASons = (4 9) - ) - (E P - E n ) = \e 2 + AM . (5) 

6 TT 

This anomaly arises, not because of an error in the pro- 
ton's energy but, because there is a nonzero Coulomb 
contribution for a (dissociated) neutron. In principle we 
should add to the above equation the contribution from 
the up-down quark mass difference. However, in the 
high-density limit, all contributions from Am are sup- 
pressed by the large Fermi momentum. Indeed, the dif- 
ference in the Fermi energy of free down and up quarks is: 
y/k F + m 2 d -^/k 2 F + m| ~(m 2 - m 2 )/2fc f . Thus, in the 
limit of very high density the total anomaly — including 
contributions from Am — becomes dominated by Eq. (|^). 
Moreover, the original mass difference between the neu- 
tron and proton (AM) "disappears" at high density be- 
cause the contributions from Am are suppressed and the 
Coulomb self-energies of the neutron and the proton are 
no longer relevant, as the quarks have rearranged them- 
selves into a uniform free Fermi gas. 

In summary, we expect that at high density there will 
be an ONS anomaly with a magnitude that is two-thirds 
that of the proton Coulomb exchange energy. Further- 
more, EM effects dominate over the contribution from 
Am and the sign of the anomaly is the same as that 
observed at normal density. 

Our earlier discussion suggests that the Coulomb en- 
ergy of pure neutron matter is nonzero. Below we focus 
on neutron matter because of the simple expectation that 
for point neutrons the Coulomb energy is zero. This may 
provide a signature of substructure. 

Since the above statements are only strictly true in 
the limit of very high density, we investigate their im- 
plications at normal density by performing a model cal- 
culation of neutron matter composed of valence quarks. 
While a model is necessary, our philosophy is to use a 
"minimal" one by demanding the following general fea- 
tures that any realistic model must posses. We require 
the many-quark wave function to (1) be explicitly anti- 
symmetric even for the exchange of quarks from differ- 
ent nucleons and (2) have cluster separability: the quark 
wave function of a nucleon removed to infinity must re- 
duce to that of a free nucleon, without any unphysical 
long-range interactions. Finally, we demand that (3) 
quarks be confined and (4) for the wave function to re- 
duce to free nucleons at low density and to a quark Fermi 
gas at high density. Perhaps, any model satisfying these 
general features can be used. 

Conventional quark potential models with two-body 
confining interactions do not satisfy cluster separability 
as they generate unphysical long-range van der Waals 
interactions between nucleons. String- flip models on 
the other hand, do satisfy the four properties described 
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above |11|,|12|. Unfortunately, we are not aware of any 
other models which both satisfy these properties and al- 
low a simple calculation. Thus, we employ the three- 
quark string- flip model discussed in Ref [fl3j. The model 
has nonrelativistic constituent quarks of mass m c of fixed 
red, green, and blue colors. A system of A nucleons is 
modeled with N = 3A quarks interacting via the follow- 
ing many-body potential: V = Vrc + Vgb + Vbr, where 
each term represents the optimal pairing of quarks. For 
example, the "red-green" component of the potential is 
defined as 



Here r^' is the coordinate of the jth red quark and Tp 
is its green partner in the neutron. The minimum is over 
all A\ permutations Pj of the set of A green quarks. A 
harmonic string potential v(r) = kr 2 /2 is used to confine 
the quarks and the Hamiltonian for the model becomes 



(6) 



(G) 



p2 

H=y A. 



V = - 



N v 2 

^— ' 2m c 



V 



(7) 



Each red quark is connected by harmonic strings to one 
and only one green and to one and only one blue quark. 
This insures that quarks will be confined into "color- 
neutral" clusters. For three quarks the model reduces 
to the well-known harmonic oscillator quark model. For 
neutron matter there is a very large number of permu- 
tations or ways to connect the strings. We employ an 
implementation of the linear sum assignment algorithm 
by Burkard and Derigs that efficiently finds the optimal 
permutation in a time proportional to N 3 jlj] . This al- 
lows Monte Carlo simulations with hundreds of quarks. 

The model has two dimension-full parameters: k and 
m c . Yet we are only interested in the harmonic-oscillator 
length b = (krric)^ 1 / 4 , as this sets the length scale for 
quark confinement. The root mean square radius of a 
nucleon is (r 2 ) 1 / 2 = 3~ 1//4 fr. Hence, to reproduce the ex- 
perimental charge radius of the proton (r 2 ) 1 / 2 = 0.86 fm 
we choose b = 1.13 fm. At the end we can rescale our 
results for other values of b. 

We are interested in simulating neutron matter. There- 
fore we assign to red and green quarks an electromagnetic 
charge of — e/3 and to blue quarks a charge of 2e/3. For 
simplicity we do not include any other intrinsic degree 
of freedom, such as spin or isospin. The electromagnetic 
self-energy of an isolated neutron is (a = e 2 = l/137) 



V° = - 



9tt (r 2 ) 1 / 2 



= -0.446 MeV 



(8) 



A simple variational wave function for the many-quark 
system has been constructed in Ref. 13 . It is given by 



* = <-m> ( -A ^ ] «j> 



(9) 



with $ a product of Slater determinants for the red, 
green, and blue quarks. In Ref. Jl3| A is a variational 
parameter characterizing the length scale for quark con- 
finement. At low density a value of A = l/v3 allows 
Eq. (||) to reproduce the gaussian wave function of a free 
nucleon. For simplicity we keep lambda fixed at A = 1/ V3 
for all densities. This insures that any change in the 
Coulomb energy of a neutron does not arise from an ar- 
tificial change in this length scale. 

We calculate the total Coulomb energy 



T/tot 

"Colli 



N 

E 



(10) 



of a system of N — 3A quarks in a box of volume V 
with antiperiodic boundary conditions. To minimize fi- 
nite size effects we use periodic distances to compute the 
quark separation. The neutron density of the system is 
p n = A/V. We use standard Metropolis Monte Carlo 
techniques to calculate the expectation value of the total 
Coulomb energy for the wave function given in Eq. (^|). 

Figure 1 shows the change in the Coulomb energy per 
neutron 



AV^j(V Cov] 



)-V° 

I y n ) 



(11) 



as a function of density for systems with N=96 and 264 
quarks. We have subtracted the neutron self-energy V® 
of Eq. (H) because this is included in the experimen- 
tal neutron-proton mass difference. We find AV to be 
nonzero. 
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FIG. 1. Change in Coulomb energy per neutron as a 
function of baryon density for pure neutron matter. The 
insert compares the model to a free Fermi gas (solid line) 
at high density. 

At normal density p n — 0.08 fm~ 3 and N = 96: 
Ay = — 78 ± 1 keV. The scale of this result suggests that 
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changes in the Coulomb energies of quarks can make a 
significant contribution to the ONS anomaly. More re- 
fined models may give results which are of the same or- 
der of magnitude, given the ratio of the nucleon size to 
interparticle spacing. Furthermore, we expect an addi- 
tional contribution from the up-down quark mass differ- 
ence Am. Our result is somewhat smaller than the total 
observed anomaly of the order 200 keV in mass 15 and 
300 keV in mass 39 pl|. Note that, for simplicity, we 
have calculated the average Coulomb energy per neutron 
rather than the self-energy of a single valence neutron. 
These quantities are expected to be similar. Indeed, in a 
free Fermi gas the average Coulomb energy per proton is 
three fourths of that of Eq. (|l|) . 

Figure 2 shows AV as a function of the nucleon root 
mean square radius or oscillator length at the fixed den- 
sity of p n — 0.08 fm~ 3 . Making the quark core of a 
nucleon smaller reduces AV, but not by much. Further, 
as the oscillator length is made very small the scale of 
the neutron self-energy V„ grows and this can increase 
AV. Of course, if the nucleon core is small one must use 
a large meson cloud to account for the full proton charge 
radius. This meson cloud, which we have not included, 
could further increase AV. 



In conclusion, we have considered a high-density limit 
of the Okamoto-Nolen-Schiffer anomaly to clarify the role 
of electromagnetic interactions (EM) and of the up-down 
quark mass difference Am. We have added a single neu- 
tron or proton to a quark gluon plasma. In this high- 
density limit we find that: (1) there is an ONS anomaly, 
(2) it is dominated by EM interactions rather than by 
Am, and (3) its magnitude is two-thirds of the pro- 
ton Coulomb exchange energy. We find an attractive 
Coulomb exchange energy for an added neutron because 
of the neutron's charged constituents. This suggests that 
the ONS anomaly could be closely related to the nucleon 
substructure. We use a minimal string- flip quark model 
to calculate the Coulomb energy of pure neutron mat- 
ter. The model wave function is fully anti-symmetric 
and satisfies cluster separability and quark confinement. 
At normal density, we find a nonzero Coulomb energy 
for neutron matter that could make a significant contri- 
bution to the ONS anomaly. 

This work was supported in part by DOE grants 
DE-FG02-87ER40365, DE-FC05-85ER250000, and DE- 
FG05-92ER40750. 
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FIG. 2. Change in the Coulomb energy per neutron as 
a function of the nucleon root mean square radius at the 
fixed neutron density of p n = 0.08 fm~ 3 . 

One should extend our results by using more elabo- 
rated quark models. It is important to study models 
with more intrinsic spin and flavor degrees of freedom 
along with more complete treatments of color. However, 
in these more complete models we still expect an ex- 
change or dynamical correlation between quarks associ- 
ated with the nucleon's hard core. This correlation could 
lead to a nonzero Coulomb energy for neutrons. Note 
that we have used harmonic oscillator confining strings. 
Thus, our wave function has gaussian tails. Linear con- 
finement may increase the tails and this should enhance 
the Coulomb exchange energy. 
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